Effective Potential in Glass Forming Liquids 
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The effective potential formalism is applied to glass forming liquids, choosing a coupling potential 
such that the "order parameter" , conjugated to the coupling strength, is the mean square displace- 
ment of the particles from their position in the quenched reference configuration. The potential is 
linear in some interval of its argument, signaling the coexistence of two phases in the system, one 
with low and one with high mean square displacement. Within this formalism, one can also compute 
the free energy of metastable states and their complexity. 

PACS numbers: 05.20.-y, 64.70.Pf 



The phase space of glassy systems is believed to break 
up, below some temperature T c called "mode-coupling 
temperature", into an extensive number of components, 
called "free energy valleys" or "metastable states" , sepa- 
rated by high free energy barriers. Lowering the tem- 
perature, the number of valleys accessible to the sys- 
tem becomes lower and lower, until at some temperature 
Tk , called "Kauzmann temperature" , it becomes non- 
extensive. The existence of metastable states in the tem- 
perature range Tk < T < T c can be evidenced by the 
study of a system coupled to a quenched configuration 
[1-5]. Consider a system described by the variables x, 
and characterized by the Hamiltonian H(x). Now switch 
on a coupling potential Tt e (xo,x) = —eQ(xo,x), where 
xq is a configuration of the system thermalized at tem- 
perature T, and then kept fixed (quenched), and Q(xo, x) 
is a measure of the "similarity" between the configura- 
tions xq and x (overlap). The free energy per particle 
/(e) of the coupled system should be independent from 
the configuration xq in the thermodynamic limit (self- 
averaging). The effective potential is defined as the Leg- 
endre transform v(q) = max e [/(e) + eq], and corresponds 
to the free energy of a system constrained at a fixed over- 
lap Q(xq,x) = Nq from the reference configuration xq, 
where N is the number of particles. 

The effective potential v(q) was computed in mean field 
p-spin models [2, 3], and in glass forming liquids using 
the hypernetted chain approximation [4]. In both cases 
it was found that, besides the thermodynamically stable 
minimum at q = Q, below temperature T c a secondary 
minimum appears, at high values of the overlap, signal- 
ing the presence of metastable states. The free energy 
difference between the two minima of the potential cor- 
responds to the logarithm of the number of free energy 
valleys, that is the complexity, that vanishes at the ther- 
modynamic glass transition Tk- Below Tk, the thermo- 
dynamic properties of the glass can still be evaluated by 
studying m coupled replicas of the system, and analyti- 
cally continuing the free energy to values m < 1 [1, 6-8]. 
A recent review can be found in Ref . [9] . 

It is important to point out that, in finite dimensional 
systems with short range interactions, free energy bar- 
riers are finite below T c due to nucleation effects (that 



is activated or hopping processes), and diverge only at 
Tk- For temperatures Tk < T < T c , the time needed to 
escape from a free energy valley is large but finite. More- 
over the coexistence of different phases, implies that the 
effective potential has to be a convex function of its ar- 
gument. This means that v(q) cannot have a secondary 
minimum, and the dynamical transition T c cannot be 
sharply defined. 

In this paper, I will apply the effective potential for- 
malism to simple liquids, choosing a particularly conve- 
nient coupling potential. The potential is such that the 
"order parameter" , conjugated to the coupling strength 
e, is the mean square displacement (m.s.d.) of the par- 
ticles from their position in the quenched reference con- 
figuration [10]. Let us consider a liquid composed by N 
monoatomic particles of mass m enclosed in a cubic box 
of side L, interacting with a potential energy U (r), where 
r is a 3./V-dimensional vector representing the positions 
of the particles, with periodic boundary conditions. The 
equilibrium free energy per particle, at temperature T, is 
given by (setting ks — 1) 



/eqC^-^ln 



d 3N r 
M\ 3N 



(1) 



where the integral is done over an hypercubic box of vol- 
ume V N , with V — L 3 , and A = h^2n/mT is the ther- 
mal wavelength. Let us now consider a (quenched) con- 
figuration ro thermalized at temperature T, and define 
the coupling potential 



n e {r ,r) = -Tin 



(e/eo )¥]r e -/^-^oi 2 



(2) 



where e > is the coupling strength, and eo = 
7rT(p/e) 2 / 3 , where p — N/V is the density of the liq- 
uid. The operator T runs over all the discrete trans- 
lations of step L in the 37V-dimensional space, and the 
operator P runs over all the permutations of the par- 
ticles. The configuration PTro runs therefore over all 
the configurations "thermodynamically identical" to Tq. 
It may be said that configuration Tq is "partially an- 
nealed", because those identical configurations are ef- 
fectively summed over in the partition function. When 
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e — > the potential H e (r ,r) vanishes, as can be seen 
doing the substitution 

P,T C_> J 

It can be seen, analogously, that ^ s -(r ,r) vanishes as 
well when e — > 0. On the other hand, when e > 0, the 
potential favors the configurations that, apart from a per- 
mutation of the particles, or a discrete translation of step 
L, are "near" to the configuration ro- We can now eval- 
uate the coupled free energy as 



d 3N r 
Nl\ 3N 



-f3[U(r)+n c (r ,r)} 



d 3N r 

^3N 



-(3U(r)-pe\r-r„\ 2 



(4) 



where in the last member U(r) has been extended by 
periodicity outside the domain < r» < L. The free 
energy /(e) can be also written as 



/(e) 



T 

'N 
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(5) 



where the effective potential v{r 2 ) is defined by 



w(r2)= -W^ (r2 -^ |r - r ° |2) 



The derivative of the effective potential with respect to 
r 2 can be expressed as 



dv{r 2 ) 3T 



dr 2 



2r 2 



[$(r 2 ) 1] , 



(7) 



where $(r 2 ) = [3{(r — r ) • V£/) r 2/3A, and the average 
(• • -) r 2 is done at fixed m.s.d. r 2 from the configuration 
ro, with the Boltzmann weight e~P u ( r \ From Eq. (6), 
being ro thermalized at temperature T, it can be derived 
that lim r 2^ [v(r 2 ) + (3T/2) In(2e7rr 2 /3A 2 )] = u cq (T), 
where u eq (T) is the mean potential energy per particle. 
The effective potential can then be computed by integra- 
tion of t9(r 2 ), 



v (r 2 ) = u cq (T)+ — 



In 
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It is reasonable to believe that /(e), v(r 2 ) and $(r 2 ) 
are self-averaging quantities, that is in the thermody- 
namic limit they do not depend on the configuration ro, 
but only on the temperature at which ro is thermalized. 
The potential v(r 2 ) can be interpreted as the free energy 
of a system constrained to be at m.s.d. r 2 from the ref- 
erence configuration ro. Actually, this is true only if r 2 



is not too large. Otherwise, many "thermodynamically 
identical" configurations, differing only by a permuta- 
tion of the particles or a discrete translation of step L, 
will be contained in the hypersphere jj\r — r | 2 = r 2 , 
and some "Gibbs correction factor" will be needed in the 
configurational integral in Eq. (6). The effective poten- 
tial v(r 2 ) and the free energy /(e) are related by the 
Legendre transform 



v(r 2 ) 



max 



/(e) + T ln(e/e ) 



er 



(9) 



Note that /(e) + In (e/eo) is a convex function. Being 



lmw,/(e) = / oq (T), and lim e ^ %(e) = [11], from 
Eq. (9) one then finds 



lim 

* 2 — >oo 



«(r 2 ) + ^ln(r 2 /r 2 nax ) 



/eq(T), 



(10) 



where = 3/(2e 1 / 3 7T/9 2 / 3 ). This equation has a sim- 
ple interpretation: when r 2 is large, the hypersphere 
^■|ro — r| 2 = r 2 contains every configuration of the liq- 
uid (modulo a permutation or a translation of step L) a 
number of times given by the ratio between the volume of 
the hypersphere and the volume V N /N\ of the "physical" 
phase space, that is (r 2 /r 2 ^) 31 */ 2 . This is the "Gibbs 
correction factor" , by which the configurational integral 
in Eq. (6) should be divided, when r 2 > r^^, to obtain 
the equilibrium free energy. Therefore, the effective po- 
tential v(r 2 ) represents the constrained free energy only 
(6) when r 2 < rj^. On the other hand, when r 2 > r^^, 
the liquid is essentially not constrained, and its free en- 
ergy is f eq (T). Putting together Eq. (8) and (10), the 
equilibrium free energy can be computed, obtaining 



/ eq (T) = u eq (T) + T 



\n{p\ 3 /e) + -h{r 2 )d\nr 



■ (11) 



In a "mean field" picture, it may happen that the effec- 
tive potential v(r 2 ) has a metastable minimum at some 
"Edwards- Anderson" m.s.d. r 2 A . As it was mentioned 
before, in a finite dimensional system with short range 
interactions, the free energy barriers remains finite when 
the number of particles goes to infinity, due to nucle- 
ation effects. Phase coexistence then implies that v(r 2 ) 
is a convex function [12]. Nevertheless, at temperatures 
near or below the mode-coupling temperature T c , if we 
let the system evolve from the configuration r , it will 
remain trapped for a very long time at some m.s.d. r 2 A 
from ro, that corresponds to the plateau in (r 2 (t)) ob- 
served in the free evolution. If the time needed to escape 
from the metastable state is longer than the time needed 
to thermalize inside it, one can effectively define the free 
energy of the metastable state as the constrained free 
energy v(r 2 A ). The complexity, defined as S = In A/", 
where Af is the number of metastable states, can then be 
computed as the difference between the constrained free 
energy v(r 2 A ), and the equilibrium one f eq (T), divided by 
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FIG. 1: Open circles: function i?(r 2 ) in LJBM for tempera- 
ture T = 0.5. The Edwards- Anderson m.s.d. is chosen (some- 
what arbitrarily) to be r| A = 0.04, that corresponds to the 
maximum of $(r 2 ), and is also close to the plateau in (r 2 (t)) 
observed in the free evolution. The interval in which points 
are well fitted by the dashed line is possibly characterized by 
phase coexistence. The critical coupling, corresponding to the 
(putative) first order transition, is e c = 1.05 ± 0.10 for this 
temperature. 



FIG. 2: Solid line: effective potential v(r 2 ) in LJBM for tem- 
perature T = 0.5. The potential is computed by Eq. (8), 
integrating a function that interpolates the measured values 
of i}(r 2 ), and setting h = m = 1. The dotted line corre- 
sponds to / q(T), computed by Eq. (11) corrected for mix- 
tures (see text). As in Fig. 1, the interval fitted by the 
dashed line is possibly characterized by phase coexistence. 
The difference [v(r EA ) — / eq (T)]/T gives the complexity, that 
is E = 0.55 ± 0.03 for this temperature. 



the temperature. We are supposing here that the "con- 
strained potential energy" (U(r)) r ^ /N coincides with 
the equilibrium one it eq (T), which is in principle assured 
only if r| A is a real minimum of v(r 2 ), so that metastable 
states have infinite lifetime. 

When temperature is low enough, a first order transi- 
tion line in the e — T plane is expected to appear. When 
this happens, there will be an interval in r 2 , correspond- 
ing to the coexistence of two different phases in the sys- 
tem, where the effective potential will be a linear function 
of r 2 , namely i>(r 2 ) = vo — e c r 2 , where e c is the critical 
value of e for that temperature. In terms of $(r 2 ), in the 
coexistence interval it will be d{r 2 ) = 1 — r 2 . 

This formalism can be applied to a hard sphere liquid 
as well. In that case, evaluating ((r — r ) • VU) r 2 as a 
temporal average, and being the time integral of the force 
— VC/ equal to the momentum variation Ap, the function 
$(r 2 ) can be computed as 



t?(r 2 ) 



lim 



ZNAt 4^ 



(r -n)- Ap, 



(12) 



where the sum is over all the collisions having place in the 
interval of time At, Ti and Api are the coordinates and 
momentum variations at the time of the i-th collision, 
and the liquid is constrained to move at fixed m.s.d. r 2 
from the configuration r . 

I have computed the function i9(r 2 ) in a well known bi- 
nary mixture of particles interacting by a Lennard-Jones 



potential (LJBM) [13], namely a mixture of 80% A par- 
ticles and 20% B particles, with pair potential v a p(r) = 
4:e al3 {(<T a [i/r) 12 - {a al3 /r)% where e AA = 1, &AA = 1, 
£ab — 1-5, <jab = 0.8, ebb = 0.5, and <jbb = 0.88. I 
simulate the system at density p — (0.94)~ 3 , where the 
mode-coupling temperature is T c = 0.435 [13]. In the 
case of mixtures, in previous equations one has to substi- 
tute N\ with NaINb*-- Moreover e has to be multiplied 
by the factor (x^Xg B ) 2 ^ 3 , where xa and xb are the frac- 
tions of A and B particles, and has to be divided 
by it. Therefore for this system one has r^ ax = 0.4220. 
Finally, the free energy in Eq. (11) gets an extra contri- 
bution T(x A In xa + xb In ig ) . 

The number of particles simulated is N = 1000, in a 
box of side L = 9.4, with pair potential truncated at a 
distance 2.5a a fj and shifted. First of all, a configura- 
tion ro is thcrmalized at temperature T, and kept fixed. 
Then, the system is cloned, copying configuration rrj into 
new variables r. Subsequently, configuration r is allowed 
to evolve, until the m.s.d. j^\r — r \ 2 is equal to r 2 . At 
this point, while r is still kept fixed, r is simulated con- 
straining the m.s.d. to remain constant, by means of the 
RATTLE algorithm [14]. It is then necessary to "ther- 
malize" configuration r, by letting it evolve onto the hy- 
persphere for a sufficient time. Finally, a time average of 
(3[(r - r ) • VU]/3N is computed to evaluate tf(r 2 ). This 
procedure can be repeated many times, with different 
and independent configurations ro, to average "over the 
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FIG. 3: Open circles: distribution of the square displace- 
ments of the particles, when the m.s.d. is constrained to be 
(r 2 ) = 0.25. Solid lines: two Gaussian distributions, with 
m.s.d. 0.053 and 0.73, and weights 0.61 and 0.23. 

disorder" , and evaluate the errors as standard deviations. 

In Fig. 1, the function d(r 2 ) for temperature T = 0.5 
is plotted. It shows a maximum at r 2 = 0.04, that is 
very close to the plateau observed in (r 2 (t)), when the 
system is allowed to evolve freely. This can be then iden- 
tified with the Edwards- Anderson m.s.d. r 2 A . The inter- 
val where #(r 2 ) is well fitted by the function 1 — r 2 
corresponds to the coexistence of two phases in the sys- 
tem. The measured values of $(r 2 ) were interpolated by 
a smooth function, that was then integrated to obtain 
the effective potential v(r 2 ) by Eq. (8), which is plotted 
in Fig. 2. Here, the interval of phase coexistence is sig- 
naled by the fit with the function u — e c f 2 • Note that the 



m.s.d. marks the crossover between the constrained 
liquid for r 2 < r^^, whose free energy is v(r 2 ), and the 
"unconstrained" liquid for r 2 > r^^, whose free energy 
is / eq (T). The complexity for this temperature is equal 
to S = 0.55 ± 0.03, a value compatible with the one 
measured in the inherent structure formalism within the 
"harmonic approximation" [15]. 

When the system is homogeneous, the distribution of 
the square displacements of the particles is expected to 
be Gaussian, P(r 2 ) oc r 2 cxp(— 3r 2 /2(r 2 )) (here we de- 
note by r 2 the square displacement of a single particle, 
and by (r 2 ) the m.s.d.) On the other hand, when two 
phases coexist in the system, we expect that the distri- 
bution is given by the superposition of two Gaussians, 
plus a contribution from the interface. This is shown 
in Fig. 3, where the distribution measured in the sys- 
tem constrained at m.s.d. (r 2 ) = 0.25 is plotted. On 
the average, a fraction 0.61 of the particles belong to the 
"bound phase", with m.s.d. 0.053, and a fraction 0.23 to 
the "unbound phase", with m.s.d. 0.73. The remaining 
fraction of particles makes up the interface between the 
two phases. This "non-Gaussian heterogeneity" shows 
up also in the free evolution of the system [16]. 

In conclusion, I have applied the effective potential for- 
malism to glass forming liquids, revealing the presence of 
a first order transition in the e — T plane, between two 
phases characterized by low and high mean square dis- 
placement with respect to a quenched reference configu- 
ration. It would be interesting to understand how much 
of this picture is relevant to the dynamical behavior of 
the system. Indeed, near or below the mode-coupling 
temperature T Cl the relaxation times may be connected 
to the free energy barrier that must be overcome to nucle- 
ate a "unbound bubble" inside the sea of bound particles. 
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